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Abstract 



Massive AGN-driven outflows are invoked by AGN-galaxy co-evolutionary models to suppress both star formation and black hole 
accretion. Massive molecular outflows have been discovered in some AGN hosts. However, the physical properties and structures of 
these AGN-driven molecular outflows are still poorly constrained. Here we present new IRAM PdBI observations of Mrk 231, the 
closest quasar known, targeting both the CO(l-O) and CO(2-l) transitions. We detect broad wings in both transitions, which trace a 
massive molecular outflow moving with velocities of up to 800 km/s. The wings are spatially resolved at high significance levels (5- 
1 lcr), indicating that the molecular outflow extends to the kpc scale. The CO(2-1)/CO(1-0) ratio of the red broad wings is consistent 
with the ratio observed in the narrow core, while the blue broad wing is less excited than the core. The latter result suggests that 
quasar-driven outflow models invoking shocks (which would predict higher gas excitation) are inappropriate for describing the bulk 
of the outflow in Mrk 231. However, we note that within the central 700 pc the CO(2-1)/CO(1-0) ratio of the red wing is slightly, but 
significantly, higher than in the line core, suggesting that shocks may play a role in the central region. We also find that the average 
size of the outflow anticorrelates with the critical density of the transition used as a wind tracer. This indicates that, although diffuse 
and dense clumps coexist in the outflowing gas, dense outflowing clouds have shorter lifetimes and that they evaporate into the diffuse 
component along the outflow or, more simply, that diffuse clouds are more efficiently accelerated to larger distances by radiation 
pressure. 
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1. Introduction 



The ubiquitous discovery of "relic" supermassive black holes 
(SMBH) in local bulges and of the tight correlations be- 
tween their masses and bulge properties (e.g. Ferrarese & Fordl 
I2005L iMarconi & Huiii l2003t iGultekin et alJ [2009t and refer- 
ences therein) indicates that nuclear activity likely plays an 
important role in galaxy evolution. Considerable attention has 
been drawn to powerful AGN-driven outflows as they may pro- 
vide a very efficient mechanism enabling the nuclear activ- 
ity occurring on sub-parsec scales to influence the surround- 
ing host galaxy. The effect of these AGN winds on the in- 
terstellar matter (ISM) in the galaxy is often referred to as 
"AGN feedback". More specifically, AGN feedback is invoked 
by most current galaxy evolutionary models to explain the 
shortage of very massive galaxies in the local Universe, their 
red colors, and low gas content, as well as the Mbh ~ cr re- 
lation (e.g. Silk &Rees] 1 19981 



Fabian 1999, Kins & Pounds 
2003: lFabian etal.ll2 006 Bald rv et al.ll2004lGranato et alJ 2004: 
Di Matteo et alJl2005t lElvisll2006t lLapi et alJl2005b lKing( [2005: 



Galaxies: individual: Mrk 231 - quasars: general - Radio lines:ISM - ISM: 



Menci et alJ [20061 120081: ICroton et all l2006t iNaravanan etail 
2008HHopkins & Elvisll2010l) . 



There are two possible "flavors " of AGN fee dback: feedback 
associated with the "radio-mode" (ICroton et al.l l2006) and feed- 
back associated with the "quasar-mode". In this work, we fo- 
cus only on the second one, which is the impulsive feedback 
exerted by AGNs during their bright active phases. One of the 
most common scenarios in theoretical models of quasar-mode 
feedback is that nuclear fast winds, generated by the AGN radi- 
ation pressure, deposit the energy produced by the nuclear ac- 
tivity into the ISM, compressing the ISM into a shock wave 
("blast wave") that eventually blows the cold gas reservoir, avail- 
able for both star formation and SMBH accretion, out of the 



host galaxy (lLapi et alJ 12005b iMenci et alJ I2006L l2008t iKm" 
120101 IZubovas & KjngH 2012). Models generally assume that 



these o utflows also contain a significantly large molec ular com- 
ponent (IZubovas & Kindl2012llNaravanan et ap2008h . If AGN- 
driven outflows can really expel large quantities of molecular gas 
(which constitutes the bulk of the ISM in the central regions and 
is the medium out of which stars form), then this would be an 
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extremely effective mechanism for quenching star formation in 
massive galaxies. 

Major observational breakthroughs in this field have been 
achieved by revealing massive and energetic quasar-driven 
outfl ows in several galaxies, both locally and at high red- 
shift dFischer et al.l201(A | Feruglio et al.l2010tlRupke & Veilieuxl 
201 It lAalto et al.1 2012t Tsturm et alJl201 it lAlatalo et alJl201 it 
Greene et al.ll20 1 It iNesvadba et al.ll20 1011201 It lAlexander et al ' 
20 10: iFarrah et alJl2012t ICano-Dfaz et alJl2012l iMaiolino et al 



2012) (see also the review by Fabian 2012). Within this con- 
text, Mrk 231, the closest quasar known, is the object in which 
a massive quasar-driven outflow was first discovered and that 
has been studied in the greatest detail ([Fischer et aL 2010t 



iFeruglio et al]l20Tot lRupke & Vei lleuxll201 It lAalto et alJl2012h 



Mrk 231 provides a unique, nearby laboratory for investigat- 
ing AGN feedback, showing all of the typical features expected 
for a quasar transiting from the obscured, merger-driven and 
IR-luminous accretion phase, accompanied by strong circum- 
nuclear star for mation, to the un obscured, standard phase of 
a 'blue' quasar (iLipariet all 1*2009 and references therein), by 
expelling obscuring gas and dust. It is indeed the most lu- 
minous ul tra-luminous infrare d galaxy (ULIRG) in the local 
Universe (Sa nders et al.l 120031) . with infrared (IR) luminosity 
L/s(8-1000pm) = 1.33xl0 46 erg s _1 , and also hosts a very pow- 
erful (Lh„i. agn = 1.1 x 10 46 erg s _1 ; e.g.. iRupke & Veilleuxl 
1201 lb . low-ionization broad absorption line (LoBAL) quasar 
(i.e. whe re the blue-shifted UV line s span velocities up to ~ 8000 
km s'^. lBoroson & Meversl(ll992l) argued that LoBAL quasars 
are young, heavily enshrouded AGNs, where the cocoon of gas 
an d dust has a l arge c overing factor, as supported by the results 
of iBraito et al.l (120041) . who found that Mrk 231 is a very ob- 
scured X-ray source with a column density of Nh ~ 10 24 cm" 2 . 
Mrk 23 1 also displays a distorted optical morphology, which im- 
plies that its host galaxy i s in an advanced stage of a merging 
process (Lipa ri et al.ll2009l) . The circumnuclear starburst in this 
source is very young (< 120 My r), with an estimated star forma- 
tion r ate of SFR * 200 M G yr 1 dTavlor et al.lll999tlDavies et al l 
l200l . 

A Herschel PACS spectrum of Mrk 23 1 has revealed a re- 
markably massive molecular outflow traced by prominent P- 
Cygni OH pro files at 65, 79, and 1 19 fim with velocity shi fts of 
~ 1000 km s" 1 dFischer et alJ2010h . In lFeruglio et alJlEOlOl) . we 
reported the detection of the molecular outflow based on IRAM 
PdBI observations, which revealed the broad (FWZI~1500 km 
s" 1 ) wings of the CO(1-0) transition, that are marginally resolved 
with an extension of about ~ 1 kpc, yielding an estimated mass- 
outflow rate of ~700 M yr" 1 , i.e., far higher than the ongoing 
SFR in Mrk 23 1 . On the basis of kinetic power arguments, we 
also suggested that the molecular outflow is primarily driven by 
the quasar radiation pressure. 

These findings were later confirmed by the analysis of the 
broad wings of the Ha and Na I D absorption line via IFU op- 
tical spectroscopy, which spatially res olved the outflow on th e 
same scales a s traced by the CO wings (Rup ke& Veilleuxl2011l) . 
Furthermore. lAalto et al.l(l2012l) also detected broad wings in the 
HCN(l-O), HCO + (1-0), and HNC(l-O) emission lines of Mrk 
231, which are tracers of the high-density molecular gas, and 
confirmed the extent of the HCN(l-O) broad emission out to 0.7 
kpc. They suggested that the molecular outflow in Mrk 231 is 
clumpy and dominated by the dense phase, and has an enhanced 
HCN abundance, which may indicate either that it contains a 
shocked medium in the outflow and/or a chemistry influenced 
by the AGN. Although these characteristics ensure that Mrk 231 
is an exception in the local Universe, they are expected to be 



quite common at high redshift, where Mrk 231 -like objects are 
thought to play a major role in the forma tion of present-day , 
red massive quiesce nt elliptical galaxies (Hop kins et al.l 12008: 
ICattaneo et al.l 2009). We note, in particular, that massive molec- 
ular outflows have been revealed in additional local ULIRGs 
through the detection of prominent OH P-C ygni profiles in their 
far-IR Herschel spectra (Stu rm et aD201 lb . 

In this paper we present the analysis of a new set of IRAM 
PdBI observations of the CO(1-0) and CO(2-l) transitions in 
Mrk 23 1 . The aim of this study is to improve the previous results 
achieved by Feruglio et al. (2010), by spatially resolving the 
molecular outflow with high significance (5-1 lcr). Moreover, 
we combine the CO(1-0) and CO(2-l) emission line data to in- 
vestigate the physical conditions of the molecular outflowing gas 
in this unique source. 

A #o=70.4 km s _1 Mpc" 1 , Q M = 0.27, Q A = 0.73 cosmol- 
ogy is adopted throughout this work. 

2. Observations 

Mrk 231 was observed in CO(1-0) and CO(2-l) with the IRAM 
millimeter-wave interferometer on Plateau de Bure (IRAM 
PdBI) between June 2009 and November 2010. The CO(2-l) 
observations presented allow us to investigate the broad wings 
in Mrk231 for the first time in this transition. The new CO(l- 
0) instead allows us to improve on the previous Feruglio et al. 
(2010) observation by significantly increasing the signal on long 
baselines and, therefore, allowing us to investigate in more de- 
tail the extension of the outflow in this transition. Table Q] pro- 
vides a technical description of these observations: the CO tran- 
sition and the corresponding redshifted frequency (z =0.042 17), 
the dates of the observations, the array configurations and the 
number of antennas used, the on-source integration times, the 
beam sizes, and the relevant references. 

The data were calibrated by using the CLIC software of the 
GILDAS package Q. The main flux calibrator at 110.697 GHz 
was MWC 349, whose flux at this frequency is 1 .27 Jy; the prin- 
cipal flux calibrators at 221.210 GHz were instead 0923+392 
and 3C84, whose fluxes at this frequency are 3.11 Jy and 7.71 
Jy, respectively. The absolute flux calibration obtained is typi- 
cally precise t o better than 10% at 110.697 GHz and ~ 20% at 
221.210 GHz dCastro-Carrizo & Nerill2010h . After calibration, 
old (2009) and new (2010) data of the CO(1-0) emission line 
were merged together. 

Both the narrow-band and the wide-band (WideX) correla- 
tors offered by the PdBI were exploited in our measurements. 
The WideX provides a spectral resolution of 1.95 MHz over 
its full bandwidth of 3.6 GHz and is available in parallel to 
the narrow -band correlator. The narrow -band correlator, w hose 
maximum signal bandwidth is 2 GHz, was configured to have a 
spectral resolution of 2.5 MHz. We rebinned the uv-tables ob- 
tained with the narrow -band correlator to give frequency chan- 
nels of 12.5 MHz for CO(1-0) and 25 MHz for CO(2-l). The 
wide-band data, thanks to the larger number of line-free channels 
available, were exploited to produce uv-tables of the continuum 
emission (sampled from channels with velocities |v| > 1200 km 
s" 1 ), which was then subtracted from both the CO(1-0) and 
CO(2-l) narrow-band and Widex data. 

Data imaging, cleaning, and analysis were performed using 
the MAPPING software, which is included in the GILDAS pack- 
age. The field of view of the maps is 2.7x2.7 arcmin 2 at 1 10.697 
GHz and 1.3x1.3 arcmin 2 at 221.210 GHz. We extracted the 



http://www.iram.fr/IRAMFR/GILDAS 
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Table 1: Description of the IRAM PdBI CO observations Table 2: Integrated fluxes 



Line 


Date(s) 


Conf. 


On source 


Beam 


Ref. 


(freq.) 




(no. of ant.) 


time 


(arcsec) 




CO(l-O) 


June-Nov. 


C+D 


20hrs 


3.2x2.8 


Feruglio 


(110.607 GHz) 


09 


(5 ant.) 






et al. 2010 




Oct. 10 


C 


7.4 hrs 




This work 






(5/6 ant.) 








CO(2-l) 


Sep. -Nov. 


C+D 


4.2 hrs 


1.6 x 1.3 


This work 


(221.210 GHz) 


10 


(5/6 ant.) 









C0(l-0) CO(2-l) 




-1009-500 500 1000 -1000-500 500 1000 
velocity (km/s) velocity (km/s) 

Figure 1 : Continuum-subtracted IRAM PdBI spectra of the 
CO(1-0) and CO(2-l) emission lines. These spectra are extracted 
from circular apertures with diameters of 6 arcseconds. Top pan- 
els: full flux scale. Bottom panels: expanded flux scale to high- 
light the broad components. The two emission lines were simul- 
taneously fitted with four Gaussian functions (red profiles, see 
text) tied to have the same centers and widths for the two transi- 
tions, to reproduce the three main components of each line (i.e. 
core, blue, and red broad components). The blue contours repre- 
sent the results of the fit. 



CO(1-0) and CO(2-l) spectra from the cleaned and continuum- 
subtracted narrow-band data cubes with an aperture of a diame- 
ter of 6 arcseconds. The noise levels per channel of the CO(1-0) 
and CO(2-l) spectra are 0.4 mJy and 1.2 mJy, respectively. We 
verified that the noise distribution is Gaussian by analyzing the 
distribution histogram of the fluxes. 

3. Results 

Figure [T] shows the IRAM PdBI spectra of the CO(1-0) and 
CO(2-l) emission lines extracted from an aperture of 6 arc- 
sec. The broad com ponent, discovered in the (1-0) transition by 
Feru glio et al.l (l2010h is also clearly present in the (2-1) tran- 
sition. The spectra in Figure Q] are fitted by using two narrow 
Gaussians to fit the core of the CO lines and two broad Gaussians 
to fit the wings. In the fitting, the centers and widths of the com- 
ponents are constrained to be identical for the two CO transi- 
tions, but their relative strengths are allowed to vary. The need to 



Fluxes from uvfits in velocity intervals 





Velocity range 
(km s" 1 ) 


F[CO(1-0)] 
(Jykms- 1 ) 


F[CO(2-l)] 
(Jy km s" 1 ) 


Nairow core 
Blue wing 
Red wing 


(-300 
(-1000 
(400 ■=■ 


f 300) 
-400) 
1000) 


(64.39 ±0.13) 
(1.84 ±0.16) 
(2.32 ±0.17) 


(194.00 ±0.50) 
(3.90 ± 0.38) 
(7.08 ± 0.84) 


Total flux 


(-1000 


r 1000) 


(75.98 ± 0.22) 


(241.32 ±0.66) 


Fluxes of Gaussian components 


Component 


Velocity 
(km s~>) 


FWHM 

(km s~') 


F[CO(1-0)] 

(Jykms-') 


F[CO(2-l)] 
(Tykms" 1 ) 


NaiTow core^ 

Bluer comp. 
Redder comp. 


(6.5 ± 2.0) 
(80.9 ± 1.3) 
(-42 ± 16) 
(527 ± 30) 


(179.7 ±2.6) 
(86.3 ±5.3) 
(548 ± 35) 
(276 ±71) 


(54.6 ± 1.1) 

(24.9 ± 2.5) 
(1.96 ±0.43) 


(173.2 ±3.3) 

(52.7 ± 6.6) 
(6.1 ± 1.4) 




Total flux 




(75.7 ± 1.9) 


(240.7 ±5.6) 



Notes: Note that a flux calibration systematic error equal to the 10% of 
the flux value should be added in quadrature to the statistical error 
reported with each flux measurement in this table. 
T The FWHM and the velocity of both the Gaussians fitting the narrow 
core are reported. 



use two narrow Gaussians to fit the core probably reflects that the 
gas distribution in the central rotating disk cannot be described 
through a simple, single Gaussian. 

The fluxes extracted from the cleaned cubes may be affected 
by various uncertainties associated with the cleaning process and 
aperture effects. Therefore, to estimate rigorously the fluxes of 
the components, we used the fluxes inferred directly from the 
visibility data, which are unaffected by any cleaning problem. 
More specifically, we measured the flux inferred from the uv 
amplitudes at the 25m baseline, for which the beam width is 
14 arcsec at 230 GHz, and is even larger at lower frequencies: 
this therefore ensures that all of the flux in the wings (which 
are much less extended, as discussed below) is accounted for. 
To obtain this information, we extracted, for both the CO tran- 
sitions, the amplitude uv diagrams of each spectral line com- 
ponent separately, by averaging the frequency channels corre- 
sponding to -300 <v(km s~')< 300 (line core), -1000 <v(km 
s -i)< -400 (blue wing), and 400 <v(km &~ l )<, 1000 (red wing). 
All three of the components are continuum-subtracted. In addi- 
tion, for the visibilities of the narrow core, we also subtracted the 
contribution of the "near" wings estimated within -600 <v(km 
s-!)< -400 and 400 <v(km s ')< 600. In Table|2] we report the 
fluxes of the CO lines, integrated within these velocity intervals 
in the uv-amplitude data, along with the total line fluxes esti- 
mated within -1000 <v(km s~')< 1000, by applying the same 
methocQ. The wings are detected at a confidence level ranging 
from 8<x to 14<x, depending on the transition. 

The full width at half maximums (FWHMs), velocities, and 
fluxes of the Gaussian components resulting from the fit shown 
in Figure Q] are also reported in Table |2] The fluxes of the 
Gaussian components were corrected for aperture and cleaning 
effects by rescaling their fluxes to the values obtained from the 
uv-amplitudes in the following way: we integrated the spectra 

2 Note that the errors reported in Table 2 are statistical errors only. 
These are the ones required to infer the statistical significance and com- 
pare the relative intensity of the various components within the same 
transition. However, when comparing with other observations or among 
different transitions, a systematic error of 10% in the absolute calibra- 
tion should be added in quadrature. 



3 



C. Cicone et al.: The physics and the structure of the quasar-driven outflow in Mrk 231 



shown in Figure[T]in the same velocity intervals as those used for 
the uv-amplitudes of the narrow core and the broad wings, and 
then estimated the three different scaling factors from the ratios 
of the fluxes resulting from the uv-amplitudes to the fluxes ob- 
tained by integrating the spectra of the individual Gaussian com- 
ponents. We note that the scaling factors obtained for the broad 
wings are slightly higher than those obtained for the narrow core, 
since the spectra extracted from an aperture with a diameter of 
only 6 arcsec miss part of the broad and faint components of the 
two CO emission lines. 





Figure 2: IRAM PdBI maps of the CO( 1 -0) blue (left panel) and 
red (right panel) broad wings. The size of the map is 20 X 20 
arcseconds. Contours correspond to 3<x (lcr = 0.14 mJy beam" 1 ). 
The synthesized beam size is shown at the bottom of the maps. 
The cross indicates the peak of the radio (VLBI) emission. Note 
that the peaks of the wing maps appear to be slightly offset from 
the peak of the radio VLBI but are well within the beam. 



CO(I-O) blue wing 



5 i 




50 100 150 200 

UV radius (m) 

C0(2-l) red wing 




Figure 4: Amplitude of visibilities (Jy) plotted as a function of 
the uv radius for the blue and red wings of the CO(1-0) and 
CO(2-l) emission lines (black points with errorbars). The visi- 
bilities were binned using intervals of 50 m of uv radius, cov- 
ering baselines from 10 m to 200 m. The emission from the 
wings was estimated within the velocity intervals -1000 <v(km 
s -i)< -400 and 400 <v(km s ')< 1000. We recall that for a 
point (i.e. unresolved) source the visibility amplitude is a con- 
stant function of the uv radius. The red points are the results of 
the best fits to a circular Gaussian model. 





Figure 3: IRAM PdBI maps of the CO(2-l) blue (left panel) and 
red (right panel) broad wings. The size of the map is 10 x 10 
arcseconds. Contours correspond to 3<x (lcr = 0.6 mJy beam" 1 ). 
The synthesized beam size is shown at the bottom of the maps. 
The cross indicates the peak of the radio (VLBI) emission. A 
tapering with uv taper = 50 m has been applied to both the maps; 
the tapering cuts the visibilities with uv radius < 50 m, so the 
maps appear more detailed and the synthesized beam is slightly 
smaller than in the original maps. Note that the peaks of the wing 
maps appear to be slightly offset from the peak of the radio VLBI 
but are well within the beam. 



Figures [2] and [3] show the maps of the blue and red wings 
of the CO(1-0) and CO(2-l) emission lines and the half-power 
contours of the synthesized beams. The emission from both the 
broad components of the CO(1-0) and CO(2-l) lines is spa- 
tially resolved. The peaks of the red and blue wings of both the 
CO(1-0) and the CO(2-l) emission lines are not spatially offset 
from one another. This validates the hypothesis that these high- 
velocity components are due to the motion of outflowing gas, 



which moves in roughly an axisymmetric way relative to the line 
of sight, instead of a rotating molecular disk. That a large frac- 
tion of the outflow occurs along our line of sight is supported 
by the blu eshifted OH absorp tion observed in a Herschel PACS 
spectrum dFischer et al.ll2010b . 

To quantify the extension of the wings, in Figure [4] we plot 
the amplitude of the visibilities as a function of the uv radius 
for both the blue and red wings of the CO(1-0) and CO(2-l) 
emission lines. We have binned the visibilities in baseline steps 
of 50 m. The data points in Figure [4] were fitted with both a 
point (unresolved) source model and a circular Gaussian source 
model. The resulting reduced x 1 an d the associated probability 
P values of the best fits are reported in Table [3] which also gives 
the differences A^ 2 =,y 2 (point)-^ 2 (Gaussian) and their associated 
I - P. We estimated the physical extent of the outflow in terms 
of the FWHM of the best-fit circular Gaussian model fits for the 
uv data. It should be noted that since we directly fitted the uv 
data, the fit automatically takes into account the interferometric 
beam size. The resulting sizes are reported in Table [3] We note 
that while the uvfit is generally good, in some cases (at some uv 
radii), the Gaussian fit to the amplitudes deviates by about 1 .5cr 
from the observed values, this means that the spatial distribution 
of the flux is probably more complex than the simple circular 
Gaussian approximation assumed by us. 

The results in Table ^demonstrate that both the red and blue 
wings of the CO(1-0) and those of the CO(2-l) lines are spatially 
resolved at high levels of significance (much higher than previ- 
ously reported in Feruglio et al. 2010 for the CO(1-0) thanks 
to the stronger signal on long baselines even for this transition). 
More specifically, the spatial distribution of the wings is resolved 
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Mrk231 core 

Mrk231 blue broad component 
Mrk231 red broad component 

Milky Way inner disk 
M82 outer disk / 
M82 central component ' 
NGC1266 CVC component /' 
NGC1266 broad component ' 
Mrk231 (Papadopoulos et al., 2007) 




J upper level 

Figure 5: Relative strengths of both the narrow core (black) and 
the broad wings (blue and red) of Mrk 231's CO emission lines. 
The ratios of the CO integrated fluxes were calculated using 
the data reported in the upper part of Table |2] The CO spec- 
tral line energy dis tribution (SLED) of Mrk 231 obtained by 
Papadopo ulos et al.l d2007l) is shown f or comparison, alon g with 
the CO SLEDs o f the Milky W ay dFixsen etalJ fl999h . M82 
(IWeifietal.ll20"05h and NGC1266 (Alata lo et al.ll2011h . 



with a confidence ranging from 5 cr to llcr, depending on the 
transition and the wing component. These results prove that the 
outflow extends out to the kpc scale: FWHM~1.2 kpc in the case 
of the CO(l-O) transition and FWHM-0.8 kpc for the CO(2-l) 
transition. 

Summarizing, our new data not only detect for the first time 
the broad wings of the CO(2-l) transition, and resolve them spa- 
tially, but also allow us to estimate the extension of the CO( 1 -0) 
broad wings with m uch higher accuracy than previous studies 
dFeruglio et al.ll2010l) . 



4. Discussion 

4.1. The CO(2-1)/CO(1-0) outflow excitation ratio 

The physical conditions of molecular gas in galaxies, such as 
its kinetic temperature, volume density, and optical depth, can 
be derived from the CO excitation properties and, in particular, 
from the relative strengths of the CO emission lines. The phys- 
ical conditions of the molecular gas responsible for the narrow 
components of Mrk 231's CO emission line s have been widely 
investi gated. Papadopoul os et al .1 d2007h and van der Werf et al. 
(2010) found that the rotational excitation diagram of the nar- 
row core of Mrk 23 1 's CO lines is dominated, in the lower J 
levels (J U ppe,- e [1,3]), by a component of density n(H2) ~ 10 3 
crrT 3 and excitation temperature T e [55,95] K. In principle, 
a similar analysis could be undertaken for the CO wings trac- 
ing the outflow. However, two CO transitions are certainly not 
enough to investigate the excitation and physical conditions of 
the outflowing gas. As a consequence, we simply compare in the 
following the CO(1-0)/CO(2-1) ratio observed in the wings with 



the same ratio observed in the core of the line and with other 
environments. 

Figure [5] shows the CO(1-0)/CO(2-1) ratio for the two wings 
and for the line core, compared with CO SLEDs observed in 
other galaxies, mor e specifically: the Milky Way's inner disk 
(2.5 < \l\ < 32.5) dFixsen et all 1 1 999h. the outer disk and the 
central component of M82 dWeifi et alj ! 2005), both t he cor e and 
the broad component of NGC1266 dAlatalo et"ail 1201 ll). and 
the CO core SLED of Mrk 231 obtained by Papadopoulos et al. 
(12001 . It should be noted that the error bars in Figure |5]repre- 
sent only the statistical uncertainties, since here we are primar- 
ily interested in the excitation of the wings relative to the core 
(hence the absolute calibration uncertainties are unimportant). 
However, when comparing excitation ratios with other targets, 
an uncertainty of 10% should be added in quadrature. 

The CO(1-0)/CO(2-1) ratios observed in the red wing and in 
the narrow core are consistent with each other. The blue wing 
is less excited than the core of the line. Taken at face value, this 
result does not support AGN feedback models where the outflow 
is generated by a radiation-driven shock ("blast wave models"), 
which, besides accelerating th e ISM, should greatly increase its 
tempe rature and density (e.g. lLapi et all2005t IZubovas & King 
12012b . However, it is possible that a high excitation component 
of the outflow may be revealed by high-J transitions. Therefore, 
sensitive and wide band observations of higher CO transitions 
are required to further investigate this issue. 

However, it should be noted that, since the CO(2-l) wings 
are more compact than the CO(1-0) wings, the CO(2-l)/CO(l- 

0) ratio is likely to vary radially. More specifically, the CO(2- 

1) /CO(1-0) ratio should be higher in the central region than 
the outer region. Higher angular resolution maps would be re- 
quired to properly map the CO(2-1)/CO(1-0) ratio and disen- 
tangle the beam-smearing effects. However, on the basis of the 
circular Gaussian modeling of the spatial distribution of the two 
wings in the two transitions, we estimated that within the central 
0.7 kpc (FWHM of the blue wing) the CO(2-1)/CO(1-0) ratio 
of the blue wing increases to 3.5+0.8, which is however still 
consistent with the excitation observed in the core of the lines 
(Fco(2-i)/F C 0(i-o) = 2.9). The CO(2-1)/CO(1-0) flux ratio of the 
red wing increases in the same inner region to 4.0+0.8, which 
starts to be inconsistent with the excitation in the core of the 
line, although marginally. On the other hand, in the region out- 
side the central 0.7 kpc the CO(2-1)/CO(1-0) ratio drops signif- 
icantly, below 2, well below the excitation observed in the core 
of the line. Summarizing, the outflowing gas associated with the 
blue wing does not show any evidence of higher excitation rel- 
ative to the core of the line, both globally and even in the cen- 
tral region, indicating that this outflowing gas is unaffected by 
shocks. Within the central region, the gas associated with the 
red wing does show marginally higher excitation, suggesting that 
this component of the outflowing gas may be affected by some 
shocks. 

The high Mach number inferred by Ferug lio et al.l (120101) . 
based on the outflow velocity of the molecular clouds, is also 
indicative of shocks. However, we also note that high outflow 
velocities do not necessarily translate into a large Mach num- 
ber, since the medium in which the clouds are embedded may 
be outflowing with the same velocity (as appears to be the case 
from the similar velocity of the ionized gas found by Rupke & 
Veilleux 201 1), hence resulting in a smaller differential velocity 
between the molecular clouds and the surrounding gas. 

Another important implication of our result is that the CO- 
to-H2 conversion factor in the wings may not differ too much 
from the bulk of the gas in the host galaxy traced by the nar- 
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4.2. Extension of the outflow 



Mrk231 red broad component 

- Mrk231 blue broad component 

CO(l-O) 
1.5 — I 




0.5 — 
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l°gK) [cm -3 ] 

Figure 6: Size (FWHM) of the outflow plotted as a function of 
the critical density of the corresponding molecular transition 
used as a wind tracer. Our CO(l-O) and CO(2-l) IRAM PdBI 
data are repr esented by filled c ircles, while the HCN(l-O) ob- 
servations bv lAalto et al.l (1201 2|) are denoted by filled triangles. 
Critical densities were calculated for a molecular medium at 
T — 100 K. Note that the same conversion factor of 0.828 kpc 
arcsec -1 has been applied to all data in order to allow compar- 
isons. 



row CO component (i.e. a ULIRG-like conversion factor). In 
Feru glio et al.l d2010l) . we had been conservative by assuming, 
for the outflow, a conversion factor a = 0.5 M (K km s~'pc 2 )~' 
(i.e. one-tenth of the Galactic value and 40% lower than the 
value assumed for ULIRGs); this is the lowest conversion fac- 
tor found in different locatio ns of M82, including its molecu- 
lar outflow dWeifi et al.ll2001l) . The finding that the CO excita- 
tion in the outflow does not differ significantly from (or is even 
lower) than the core of the line, tracing the bulk of the gas in 
this ULIRG, suggests that the conversion factor is also similar 
(i.e. a = 0.8 M (Kkm s 'pc 2 ) -1 ). This imp lies that the outflow 
molecular mass and outflow rate given in Ferug lio et al.l (2010) 
are underestimated, and are probably higher by a factor of about 
1.4. 

We note that i n the j et-driven outflow observed in NGC 1266 
by lAlatalo et alJ (1201 ll) the excitation of the CO broad compo- 
nents tends to be higher than in the central velocity component 
(as illustrated in Fig. [5j, although not strongly. This is not unex- 
pected. In this case, the outflow is totally driven by jet-induced 
shocks, which unavoidably heat the ISM and therefore produce 
higher CO excitation. The outflow CO excitation differences be- 
tween Mrk 231 and NGC 1266 further highlight that the driving 
mechanism in the two sources are different, although both asso- 
ciated with an AGN. 

Finally, we note tha t the M rk 23 1 SLED of the core obtained 
by IPanadopou los et all d2007l) by means of single dish obser- 
vations is slightly steeper than the one obtained by us. We as- 
cribe the difference partly to calibration uncertainties and partly 
to uncertainties in subtracting the continuum and broad wings 
with previous narrow band spectra (especially in single dish data 
where the baseline instabilities may be problematic). 



As mentioned above, our new data have allowed us to constrain 
the size of the molecular outflow traced by the CO transitions 
with unprecedented accuracy, enabling analyses that have not 
been possible until now. The results presented in the previous 
section clearly reveal that the extension of the outflow traced 
by the CO(1-0) transition is significantly larger than observed 
in the CO(2-l) and HCN(l-O) transitions. Fig. |6] shows the ex- 
tension of the outflow as a function of the critical density of the 
transition used to trace the outflow, where the red and blue sym- 
bols indicate the trends relative to the red and blue wings, re- 
spectively. There is a clear trend for the size of the outflow to 
decrease as a function of the critical density, at least in the range 
10 3 < n cr < 10 4 cm (while at higher densities it is more dif- 
ficult to assess this trend owing to the large errorbars in the size 
of the HCN wings). 

A radial decreasing density of the outflowing molecular gas 
is also consistent with the Herschel measurements by Fischer 
et al. (2010) and Sturm et al. (2011), and in particular with the 
blueshifted absorption profile of the far-IR molecular transitions 
observed in the PACS spectra (OH 79/mi, OH 119//m, 18 OH 
120//, H2O 79//m), as discussed in those papers. 

This result is consistent with the prediction of 

iNaravanan et al.l J2008), who studied the response of CO 
morphology to starburst and AGN feedback-driven winds in 
galaxy mergers. According to them, the emission from higher 
CO transitions (e.g. CO(2-l)) broadly follows the distribution of 
the CO(1-0) emission, but remains slightly more compact. This 
is because, according to their model, denser clouds (preferen- 
tially traced by higher transitions) have shorter lifetimes along 
the outflow and tend to evaporate into a more diffuse molecular 
phase, traced by the CO(1-0) transition. 

The result could also possibly support the scenario proposed 
by Hopkins & Elvis (2010), where dense clouds invested by an 
even mild AGN-driven wind develop Kelvin-Helmholtz insta- 
bilities causing them to expand and accelerate more easily by 
means of radiation pressure. In the latter scenario, one would 
expect diffuse clouds to be driven to larger distances than the 
denser clouds. However, Kelvin-Helmholtz instabilities should 
also develop shocks (although much milder than those expected 
in "blast wave models") that should show up in higher gas ex- 
citation, which is not the case for the blue wing, but may be for 
the red wing. 

Simple accel eration by radiatio n pressure on the dust in 
molecular clouds (Fabian et a l.l2009l) . without invoking any sce- 
nario associated with shocks, may not only explain the mod- 
est CO excitation in the wings, but also the size difference be- 
tween the CO(1-0) and CO(2-l)+HCN emission components. 
More diffuse clouds (preferentially traced by CO(1-0)) are in- 
deed characterized by a higher ratio between the outward radia- 
tion pressure and inward gravitational pull than dense clouds, 
which are preferentially traced by CO(2-l)+HCN emission. 
More specifically, for clouds whose dust is optically thick to 
UV/optical radiation, the ratio F ra d/F grav is inversely propor- 
tional to the column density of the cloud Nh, which is higher in 
dense clouds. As a consequence, diffuse clouds can be expelled 
to larger distances. Radiation pressure on the dust in molecular 
clouds can in principle also develop shocks, but this depends on 
whether the medium in which the molecular clouds are embed- 
ded is comoving with them. 

We have discussed so far only the global size of the emitting 
regions of the two CO transitions obtained by fitting the uv data 
with a simple circular Gaussian. At our resolution, the maps do 
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Table 3: Best- fit results of the visibility versus uv radius plots. 



Line 


Velocity range 


Point 


Gaussian 


4r 2t 


FWHM 


FWHM* 




(km s _1 ) 


source model 


source model 




(arcsec) 


(kpc) 


CO(l-O) 


[-1000, -400] 




X 2 r = 1-98 


AX 2 = 6.48 


(1.52 ±0.31) 


(1.26 ±0.25) 


blue wing 




P = 6.1E-03 


P = 0.11 


(1 -P) = 0.99 






CO(l-O) 


[400, 1000] 


xl = 4.26 


Xl = 1-74 


A x 2 = 7.54 


(1.37 ±0.22) 


(1.13 ±0.18) 


red wing 




P = 5.2E-03 


P = Q.\6 


(l-/>) = 0.99 






CO(2-l) 


[-1000, -400] 


xl = 1-69 


^=0.14 


AX 2 = 4.65 


(0.80 ±0.16) 


(0.66 ±0.14) 


blue wing 




P = QA1 


P = 0.94 


(1 -P) = 0.91 






CO(2-l) 


[400, 1000] 


X 2 r = 8.94 


^=0.19 


Ax 2 = 26.24 


(1.01 ±0.09) 


(0.84 ± 0.07) 


red wing 




P = 6.5E-06 


P = 0.90 


(1 -P) ^ 1.00 







Notes: * A^ 2 = / ^ 2 (point)- / ( : ' 2 (Gaussian) * The adopted cosmology and redshift result in a spatial scale of 0.828 kpc arcsec 



not show much structure. The only relevant structure at the 3<x 
level is possibly an extension of the CO(1-0) blue wing (Fig. 2, 
left) a few arcseconds to the north. This could imply that there is 
an association between some of the molecular outflow (a minor 
part of it) and the northern region, where the radio jet affec ts the 
neutral gas outflow as inferred bv lRupke & Veilleux (201 1) from 
the blue-shifted Na I D kinematics map. lAalto et akl d2012h de- 
tected a plum extending to the north in the HCN(l-O) red-shifted 
wing map and also inferred a link betwee n this feature and 
the je t-influenced neutral wind detected by Rupke & Veilleuxl 
(I2TJ1I . However, we propose that this association is incorrect, 
since the jet-accelerated gas is detected in terms of blueshifted 
Na I D absorption, so any association should be seen in the 
molecular blue wings, not the red wings. The possible inter- 
action between part of the molecular outflow and the radio-jet 
should however be investigated using higher-resolution data. 

Finally, we note that the size (~2 arcsec) of the (blueshifted) 
"nuclear wind", which is the region identified by Rupke & 
Veilleux (2011) where the outflow is dominated by a radiation- 
pressure drive wind, is close to the size (1.5 ± 0.3 arcsec) of the 
molecular outflow determined by us from the blue wing of the 
CO(1-0) component. 



5. Conclusions 

We have presented new broad-band IRAM PdBI observations of 
the CO(2-l) transition in Mrk 231, which is the ULIRG host- 
ing the closest quasar known. For the first time, our observations 
reveal the broad wings of the (2-1) transition extending to ve- 
locities of FWZI~1500 km s _1 , which trace the same molecular 
wind revealed by previous CO( 1 -0) PdBI observations and PACS 
far-IR spectroscopy of OH transitions. The CO(2-l) wings are 
spatially resolved, with an extension of about 0.8 kpc. 

We have also obtained new CO(1-0) observations, by sig- 
nificantly improving the signal at high angular resolution rela- 
tive to previous works, which allows us to determine more ac- 
curately the size of the molecular outflow traced by this tran- 
sition. The CO(1-0) broad wings are spatially resolved with a 
significance that is much higher (5-7 cr) than in previous obser- 
vations. We measure an extension of the CO(1-0) wings of about 
1 .2 kpc, which is significantly larger than observed for the CO(2- 
1) wings. The extension of the CO(1-0) wings is also signifi- 
cantly larger than the extension of the HCN(l-O) wings, tracing 
high density gas in the outflow detected by previous observa- 
tions. More specifically, we show that the size of the outflow 



anticorrelates with the critical density of the transitions used to 
trace the outflow. 

These results are consistent with the scenario where denser 
clouds (traced by the CO(2-l) and HCN(l-O) broad wings) have 
shorter lifetimes along the outflow and evaporate into a more 
diffuse molecular component (traced by CO(1-0)) at larger radii, 
as predicted by some feedback models. 

We also find that the CO excitation in the outflow, as traced 
by the CO(2-1)/CO(1-0) ratio of the broad wings, does not differ 
significantly from that of the gas in the bulk of the galaxy, as 
traced by the core of the CO lines. In the blue wing, the exci- 
tation is lower than in the core of the line. Taken at face value, 
this result in inconsistent with models where the molecular out- 
flow is driven by a shock wave generated by the interaction of a 
radiation pressure-driven nuclear wind with the ISM of the host 
galaxy. We favor a scenario where gas clouds are directly ac- 
celerated by the radiation pressure on dust. This scenario does 
not need to invoke shocks, hence the CO excitation is unaf- 
fected. Moreover, in this scenario low density clouds (traced by 
CO(1-0)) are accelerated more efficiently, hence reaching larger 
distances than dense clouds (traced by CO(2-l) and HCN). 
However, we also note that in the inner region (R < 0.3 kpc), 
the CO(2-1)/CO(1-0) ratio is indeed slightly higher, which is in- 
dicative of some shock contribution. 

We marginally detect an extension of the CO(1-0) 
blueshifted wing to the north, where previous studies have found 
evidence that the atomic neutral outflow, traced by blueshifted 
Na I D absorption, is influenced by the radio jet. This suggests 
that the radio jet may also contribute to the acceleration of part 
of the molecular outflow, although far less significant than the 
global molecular outflow. 
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